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ABSTRACT 


Vorticity  forms  the  basis  for  a  new  approach  to  the  anal- 
ysis of  inlet  distortion  in  axial  flow  fans  and  compressors. 
Farmer  first  formulated  this  approach  and  this  paper  repre- 
sents a  test  of  the  usefulness  of  his  theory.   A  summary  of 
recent  developments  in  the  field  along  with  a  discussion  of 
the  effects  of  distortion  on  stall  margin  forms  a  background 
for  this  complex  problem.   A  computer  program  calculates  the 
vorticity  pattern  at  the  compressor  face  using  data  that  are 
read  from  a  magnetic  tape.   The  data,  which  were  provided  by 
NASA  Lewis  Research  Center  in  Cleveland,  consist  of  eleven 
stall  events.   The  output  of  the  program  contains  tables  of 
stagnation  pressure  and  the  partial  derivatives  of  pressure 
in  the  R  and  theta  directions  as  well  as  three  maps.   The 
three  maps  are  of  pressure  and  radial  and  circumferential 
vorticity.   The  results  of  the  program  show  a  correlation  be- 
tween a  ring  of  large  positive  circumferential  vorticity  and 
stall,  leading  to  the  conclusion  that  the  stall  was  caused  by 
the  increase  in  blade  loading  induced  by  the  vorticity.   This 
conclusion  suggests  a  formulation  for  a  universal  inlet 
distortion  index  and  provides  a  basis  for  the  evaluation  of 
the  vorticity  approach  to  the  problem. 
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I.   INTRODUCTION 

The  problem  of  surge  in  axial  flow  fans  and  compressors 
in  jet  engines  is  a  serious  one.   Unsolved  it  can  result  in 
the  loss  of  a  multi-million  dollar  aircraft  and  possibly  in 
the  loss  of  the  life  of  the  pilot.   Manufacturers  are  there- 
fore forced  to  design  compressors  which  operate  with  decreased 
performance. 

During  the  early  years  of  the  use  of  axial  compressors  in 
aircraft  jet  engines,  before  high  speed  flight  was  a  common 
occurrence,  the  problem  of  surge  was  treated  as  a  steady-state 
problem  with  satisfactory  results.   However,  with  the  advent 
of  high  performance  and  high  speed  military  aircraft,  and  in 
particular  the  F-lll  which  employed  the  TF-30  turbo-fan  engine, 
[Ref.  1]  it  became  evident  that  the  problem  could  no  longer  be 
considered  to  be  steady-state.   There  were  some  attempts  made 
to  improve  the  steady-state  distortion  factor,  but  these  were 
unsuccessful.   During  the  next  few  years,  various  attempts 
were  made  to  develop  some  sort  of  a  distortion  index  using  the 
fluctuations  of  the  stagnation  pressure  as  a  starting  point. 
These  attempts  met  with  some  degree  of  success,  but  none  was 
the  complete  answer  to  the  problem. 

It  has  been  proposed  that  a  distortion  factor  which  uses 
vorticity  as  a  measure  of  instantaneous  distortion  would 
provide  insight  into  the  fluid  mechanics  of  the  problem. 
Farmer  [Ref.  2]  first  suggested  this  approach  and  developed 


approximate  equations  for  the  vorticity  at  the  compressor 
face  in  terms  of  measured  quantities,  i.e.,  stagnation 
pressure  data. 

Since  most  of  the  distortion  factors  developed  by  other 
researchers  were  the  outgrowth  of  instantaneous  maps  of  the 
pressure  distribution  over  the  compressor  face,  it  seems 
logical  to  take  that  same  route  to  the  development  of  a  dis- 
tortion index  based  on  vorticity.   Therefore  a  computer  pro- 
gram was  developed  to  calculate  the  instantaneous  vorticity 
pattern  and  produce  the  maps  of  the  pattern  on  the  printer 
component  of  the  computer.   This  paper  will  attempt  to  inter- 
pret the  output  of  this  program  and  prove  the  validity  of 
this  approach.   Suggestions  then  will  be  made  concerning  the 
formation  of  a  distortion  factor  that  will  allow  the  pre- 
diction of  surge  for  a  wide  variety  of  engines. 


II.   NATURE  OF  PROBLEM 

A.   SOME  PERTINENT  DEFINITIONS 

In  order  to  understand  the  complexity  of  the  problem,  one 
must  first  be  familiar  with  some  of  the  nomenclature  and  basic 
concepts  that  are  used  in  the  field  of  compressor  stability. 
Brimelow  [Ref.  3]  gives  an  excellent  presentation  of  the 
present  theories  on  this  subject  as  well  as  a  concise 
explanation  of  the  nomenclature. 

The  most  likely  place  for  an  engine  to  become  unstable  is 
in  the  fan  or  compressor.   These  two  components  are  said  to 
be  in  surge  when  several  or  all  of  the  blades  become  too 
heavily  loaded  and  are  forced  to  operate  in  an  aerodynamically 
stalled  condition.   If  one  examines  the  compressor  map  of  the 
J-85-GE-13  jet  engine  in  Fig.  1,  he  can  see  that  the  conditions 
for  surge  are  represented  by  a  line  on  the  compressor  map. 
This  line  is  generally  called  the  stall  or  surge  line.   The 
most  efficient  operating  condition  of  the  compressor  is  as 
close  to  this  stall  line  as  possible  without  stalling.   In 
order  to  insure  that  the  compressor  will  not  surge,  it  must  be 
run  with  a  safety  margin  between  the  operating  point  and  the 
stall  line.   The  vertical  difference  between  the  operating 
point  and  the  stall  line  is  defined  as  the  stall  or  surge 
margin.   In  today's  high  performance  aircraft  it  is  necessary 
to  have  as  small  a  surge  margin  as  possible  in  order  to  de- 
crease the  specific  fuel  consumption  and  increase  the  total 
thrust,  while  building  the  engine  to  be  as  light  as  possible. 
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FIG,  1  COMPRESSOR  MAP  J-85-GE-13 


It  appears  to  be  an  easy  task  to  define  a  small  surge  mar- 
gin and  then  run  the  engine  at  a  point  such  that  this  surge 
line  will  never  be  crossed.   It  turns  out,  however,  that  the 
compressor  map  is  only  a  representation  of  the  steady-state 
operating  conditions  of  the  compressor.   Once  that  fact  is 
known,  it  is  understandable  that  fluctuations  in  the  flow  may 
affect  the  stall  margin  and  cause  the  compressor  to  stall  even 
though  the  steady-state  operating  point  may  have  been  well 
below  the  stall  line. 

The  term  compressor  instability  may  best  be  defined  through 
visual  means.   Figure  2  is  a  hypothetical  compressor  map  with 
only  one  constant  speed  line.   On  this  line  six  operating 
points  are  labeled  PI,  P2,  Rl ,  P3,  P4 ,  and  R2.   If  the  com- 
pressor is  operating  on  the  negative  slope  of  this  line  at 
point  Rl ,  its  operation  is  stable.   It  is  stable  because  the 
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operating  point  will  tend  to  return  to  Rl  if  it  is  forced  in 
either  direction  along  the  constant  speed  line.   If  it  is 
forced  to  P2  the  mass  flow  will  be  decreased.   The  character- 
istics of  the  chamber  behind  the  compressor,  which  is  known 
as  the  receiver,  require  a  lower  pressure  ratio  at  the  re- 
duced mass  flow  rate.   The  compressor  discharge  flow  will  be 
forced  to  expand  into  the  receiver.   This  will  cause  an  in- 
crease in  mass  flow  rate  and  the  operating  point  will  return 
to  Rl.   In  a  similar  manner  if  the  operating  point  is  forced 
to  Pi,  it  will  be  forced  to  return  to  Rl .   In  this  case  the 
pressure  ratio  will  be  less  than  that  required  by  the  receiver 
at  the  increased  mass  flow  rate  and  a  back  pressure  will  re- 
sult.  The  back  pressure  will  cause  a  reduction  in  flow  rate 
and  the  operating  point  will  return  to  Rl. 
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FIG,  2  COMPRESSOR  STABILITY 
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If  the  compressor  is  operating  on  the  positive  slope  of 
the  constant  speed  line  at  point  R2 ,  the  operation  of  the 
compressor  is  unstable.   In  this  case,  if  the  operating  point 
is  disturbed  from  R2 ,  it  will  continue  to  move  in  the  direc- 
tion in  which  it  initially  was  moving  and  will  not  return  to 
R2 .   If  the  operating  conditions  move  to  point  P4,  a  decrease 
in  pressure  ratio  will  occur.   This  will  in  turn  set  up  a 
back  pressure  since  the  pressure  in  the  receiver  will  be 
greater  than  the  compressor  discharge  pressure,  and  the  flow 
rate  will  be  reduced.   The  decrease  in  flow  rate  will  con- 
tinue to  reduce  the  pressure  ratio,  and  the  operating  point 
will  continue  to  move  away  from  R2 .   If  the  operating  point 
moves  in  the  other  direction  along  the  constant  speed  line, 
the  pressure  ratio  will  increase  which,  in  turn,  will  in- 
crease the  flow  rate,  and  again  the  operating  point  will 
continue  to  move  away  from  the  design  point  R2 . 

One  can  see,  therefore,  that  the  limit  for  stable  com- 
pressor operation  is  at  the  point  where  the  slope  of  the 
constant  speed  line  on  the  compressor  map  goes  to  zero.   The 
compressor  must  operate  where  the  slope  of  the  constant  speed 
line  is  negative  and  not  be  allowed  to  reach  the  pressure 
ratio  that  corresponds  to  the  limiting  point  for  stable  com- 
pressor operation. 

With  high  performance  aircraft  the  transient  distortions 
induced  by  sources  external  to  the  engine  and  in  the  engine 
inlet  will  be  of  a  large  magnitude  and,  therefore,  more 
dangerous  to  the  stability  of  the  compressor.   The  turbulence 
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is  likely  to  be  generated  by  some  or  all  of  several  sources. 
Some  of  these  sources  are:   shock  waves,  boundary  layers,  hot 
gases  emitted  from  rockets  launched  from  the  aircraft,  and 
off -design  operation  of  the  inlet.   Since  the  stall  margin  is 
such  a  critical  variable  in  the  operation  of  a  jet  engine,  it 
is  not  feasible  merely  to  increase  the  stall  margin  to  such 
an  extent  that  time-variant  distortion  cannot  cause  the 
compressor  to  stall  without  a  disastrous  drop  in  the  overall 
efficiency  of  the  engine.   It  is  necessary  to  obtain  more 
information  on  the  source  of  the  distortion,  what  effect  the 
distortion  has  on  the  inlet-compressor  interface,  and  just 
exactly  what  type  of  distortion  is  most  likely  to  cause  a 
compressor  to  surge. 

B.   RECENT  DEVELOPMENTS  IN  THE  FIELD 

With  the  advent  of  the  new  high  performance  aircraft  it 
became  evident  that  the  method  that  was  used  for  the  evaluation 
of  the  inlet-compressor  interface  was  no  longer  effective.   At 
that  time  this  evaluation  had  been  done  wholly  by  an  analysis 
of  the  steady-state  distortion  through  the  use  of  a  steady- 
state  distortion  factor.   In  a  high  performance  aircraft  the 
inlet  flow  is  more  likely  to  be  turbulent  making  it  necessary 
to  evaluate  the  effects  of  time  variant  distortion  on  the 
compressor.   Since  a  turbulent  flow  is  a  function  of  four 
variables,  which  in  a  cylindrical  coordinate  system  are  R,  Q    , 
Z  and  time,  the  problem  became  a  much  more  complex  one.   The 
magnitude  of  the  fluctuations  in  pressure  and  velocity  in  a 
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turbulent  flow  is  such  that  they  cannot  be  considered  small 
with  respect  to  the  steady-state  values. 

In  the  light  of  turbulent  inlet  flows,  it  becomes  apparent 
why  the  steady-state  distortion  index  was  not  adequate.   It 
does  not  take  the  time  variable  into  account.   When  the  inlet 
flow  of  a  compressor  is  turbulent,  the  -steady-state  distortion 
factor  can  be  well  within  the  tolerences  for  stall,  and  yet 
the  compressor  may  stall  due  to  the  random  velocity  and 
pressure  fluctuations  in  the  flow. 

The  first  attempt  to  solve  this  problem  resulted  in  an  im- 
proved steady-state  distortion  factor.   As  might  be  expected, 
however,  that  factor  was  also  inadequate  for  the  solution  of 
the  problem.   The  amount  of  turbulence  is  often  quantified 
by  the  variation  of  root-mean-square  stagnation  pressure 
(A  Prms/P) .   Still  reluctant  to  tackle  the  turbulent  flow 
problem,  researchers  attempted  to  devise  a  method  by  which 
the  improved  steady-state  distortion  factor  could  be  used 
with  a  correlation  for  the  amount  of  turbulence  [Ref.  1], 
This  method  was  also  ineffective,  but  it  gave  some  insight 
into  the  complexity  of  the  problem. 

At  this  point,  it  was  evident  that  some  method  of  predict- 
ing the  effects  of  time  variant  distortion  must  be  developed 
to  be  used  in  the  inlet-compressor  matching  problem.   The 
method  used  was  to  introduce  a  system  of  instrumentation  that 
had  a  quick  response  and  was  able  to  measure  the  stagnation 
pressure  at  any  particular  point  in  time.   Measurements  then 
were  taken  at  small  intervals  of  time  with  probes  located  on 
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pressure  rakes  at  the  compressor  face.   The  pressures  at 
these  discrete  points  were  converted  to  a  distribution  and 
presented  in  the  form  of  an  instantaneous  pressure  map.   By 
fixing  the  variable  Z  at  the  compressor  face,  researchers 
eliminated  one  variable;  and  by  taking  instantaneous  pressure 
readings,  they  made  the  problem  two-dimensional  at  any  point 
in  time.   The  influence  of  time  was  expressed  in  the  analysis 
by  the  time  continuity  of  a  set  of  maps.   It  became  possible 
to  find  an  empirical  correlation  between  the  instantaneous 
distortion,  as  represented  by  the  instantaneous  pressure  map, 
and  surge.   In  Ref.  5  Calogeras,  Burdstadt,  and  Coltrin  found 
an  empirical  correlation  between  the  instantaneous  pressure 
distribution  and  the  onset  of  surge.   One  of  the  maps  they 
used  in  the  development  of  the  correlation  is  depicted  in 
Fig.  3.   Using  their  correlation  between  the  instantaneous 
pressure  distribution  and  surge,  they  developed  an  instan- 
taneous inlet  distortion  factor  with  which  it  was  possible  to 
predict  stall  with  some  measure  of  success.   This  factor  was 
developed  for  the  J-85-GE-13  engine  with  a  particular  inlet 
geometry  and  was  not  tested  with  either  other  engines  or 
other  inlets.   It  was  tested  at  angle  of  attack  and  proved 
quite  successful  in  that  context. 

Plourde  and  Brimlowe  [Ref.  6]  also  developed  an  instan- 
taneous inlet  distortion  factor.   This  factor  was  used  for 
the  TF-30  turbofan  engine  which  is  installed  in  the  F-lll.   It 
is  interesting  to  note  that  both  of  these  factors,  for  the  J-85 
and  the  TF-30,  arc  primarily  circumferential  distortion  factors, 
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FIG,  3  INSTANTANEOUS  PRESSURE  MAP 
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Another  interesting  point  is  that  while  each  of  these  factors 
is  useful  in  the  analysis  of  the  case  for  which  it  was  devel- 
oped, neither  would  be  useful,  except  in  a  general  way,  to 
anyone  working  on  the  inlet  distortion  problem  in  a  different 
engine. 

The  period  of  time  for  which  a  particular  distortion 
pattern  must  persist  in  order  to  cause  stall  was  found  to 
be  extremely  small.   Reference  7  states  that  for  the  TF-3  0 
this  time  may  be  approximately  three  milliseconds,  while 
Ref.  5  has  stated  that  the  time  scale  for  the  J-85  can  be  as 
short  as  one  millisecond.   Reference  8  contains  a  method  of 
calculating  the  length  of  time  for  which  a  distortion  must 
persist  in  order  to  cause  steady-state  stall.   In  general 
this  time,  depending  on  the  size  of  the  engine,  engine  speed, 
and  number  of  blades,  is  between  one  and  five  milliseconds. 
If  the  duration  of  a  distortion  is  less  than  this  time,  stall 
will  not  occur,  due  to  dynamic  suppression  of  stall,  even 
though  the  type  of  distortion  might  have  caused  a  stall  had 
it  persisted  for  a  longer  period  of  time.   Unfortunately  the 
exact  time  that  a  distortion  must  persist  cannot  be  calculated 
since  it  may  be  a  function  of  many  more  variables  than  those 
mentioned  above. 

Upon  examination  of  several  other  papers,  [9,  10,  11,  12] 
it  becomes  evident  that  the  large  amount  of  work  that  has  been 
done  on  this  problem  has  resulted  in  two  conclusions.   One  of 
these  results  is  an  empirical  relationship  between  instanta- 
neous inlet  distortion  and  stall,  while  the  other  result  is 
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lack  of  agreement  on  the  quantitative  form  of  the  relation- 
ship.  As  a  consequence  dozens  of  instantaneous  distortion 
factors  are  used  at  the  present  time.   This  lack  of  consis- 
tency can  be  realized  by  examining  Table  I  which  contains 
six  distortion  factors,  all  of  which  have  been  employed  for 
the  same  engine  at  various  points  in  time. 

TABLE  I  -  DISTORTION  FACTORS 
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Brownstein  [Ref.  13]  puts  forth  a  possible  universal  in- 
let distortion  factor  using  pressure  fluctuations  as  a  base 
for  the  index.   It  appears  to  have  solved  some  of  the  prob- 
lems encountered  by  earlier  factors  and  is  a  more  general  in- 
dex.  This  index  seems  to  be  very  promising.   Another, 
different  approach  to  the  problem  was  put  forth  by  Farmer 
[Ref.  2]  and  again  in  Ref.  8. 

C.   VORTICITY  APPROACH 

The  approach  presented  in  Ref.  2  and  Ref.  8  is  new. 

Using  Crocco ' s  theorem  it  was  possible  to  obtain  an  ex- 
pression for  vorticity  at  the  compressor  face.   Then,  with 
an  order  of  magnitude  analysis,  using  data  taken  for  the 
J-85-GE-13  jet  engine,  Farmer  was  able  to  eliminate  some  of 
the  terms  involved  in  the  original  equation  and  emerge  with 
a  simple,  easily  calculated  expression  for  vorticity  at  the 
compressor  face.   The  two  expressions,  for  radial  and 
circumferential  vorticity  are 

^      '        -1        2JP'  (1) 


UJh 


1      d_p 

*Fuzr    do  (2) 


For  a  more  detailed  explanation,  consult  Farmer's  work  in 
Ref.  1.   Due'  to  the  fact  that  vorticity  lines  do  not  end, 
i . e. , V X  UJ  =  0,  the  axial  vorticity  can  be  calculated  using 
a  control  volume.   Unfortunately  the  direction  of  the  ixial 
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vorticity  will  be  ambiguous.   As  shown  in  Fig.  4,  it  is  not 
possible  to  determine  if  the  axial  vorticity  enters  the 
control  volume  from  the  front  or  the  back. 


FIG.  '1  AMBIGUOUS  DIRECTION  OF  AXIAL  VORTICITY 


Reference  2  states  that  circumferential  vorticity  is 
likely  to  have  a  larger  effect  on  the  stall  margin  than  is 
radial  vorticity.   This  is  in  concurrence  with  other  work 
done  on  this  problem,  but  it  may  be  beneficial  to  consider 
the  effects  of  radial  vorticity  in  the  light  of  the  fact 
that  the  duration  of  a  distortion  pattern  necessary  to  cause 
stall  is  extremely  short.   It  should  be  noted,  however,  that 
radial  vorticity  in  a  rotor  is  an  unsteady  phenomenon  [Ref.  8] 
and  therefore  very  difficult  to  treat. 

This  new  approach  seems  to  have  merit,  and  its  usefulness 
should  be  tested  using  some  of  the  large  supply  of  existing 
data  that  have  been  accumulated  in  the  study  of  this  problem. 
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It  would  appear  that  some  method  of  producing  instantaneous 
vorticity  maps  similar  to  instantaneous  pressure  maps  would 
be  useful  in  the  evaluation  of  this  theory.   Also  should  the 
theory  be  confirmed  by  these  maps,  the  maps  would  be  instru- 
mental in  the  formation  of  a  distortion  index  based  on 
vorticity. 

The  high  speed  digital  computer  is  the  best  choice  for 
the  means  of  production  of  these  maps.   The  main  reason  for 
this  choice  is  the  amount  of  time  that  would  be  necessary 
for  their  production  by  hand  or  even  on  the  most  sophisti- 
cated desk  calculator  that  is  available  today.   Also,  it  is 
necessary  that  the  partial  derivatives  of  pressure  in  the 
r1  and  Q    directions  be  calculated,  and  this  can  be  done  most 
accurately  by  the  use  of  numerical  curve  fitting  routines 
that  already  have  been  written  for  the  computer. 

The  problem  then  is  to  develop  a  computer  program  that 
will  fulfill  the  following  requirements:   read-  the  data  from 
a  magnetic  tape,  define  a  function  in  terms  of  r'  and  Q    that 
will  provide  the  values  of  the  normalized  pressure  and  its 
two  derivatives  at  any  point  on  the  compressor  face,  calcu- 
late the  normalized  axial  velocity,  calculate  the  radial  and 
circumferential  vorticity  at  any  point  on  the  compressor 
face,  and  finally  print  maps  of  pressure  and  radial  and 
circumferential  vorticity  at  the  compressor  face.   These  re- 
quirements are  shown  more  clearly  by  the  flow  chart  in  Fig.  5. 

Once  the  program  has  been  developed  and  meets  the  require- 
ments mentioned  above,  the  maps  should  be  interpreted  as  to 
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the  effect  of  the  instantaneous  vorticity  distribution  on  the 
stall  margin.   It  is  hoped  that  there  will  be  some  relation- 
ship between  the  instantaneous  vorticity  pattern  and  surge 
other  than  the  purely  empirical  correlation  between  the 
instantaneous  pressure  distribution  and  surge  that  has  been 
found  before. 


MAGNETIC  TAPE 


FIG,  5  COMPUTER  PROGRAM  FLOW  DIAGRAM 


Farmer,  Iverson  and  Fuhs  [Refs.  8  and  14]  have  presented 
one  such  logical  relationship  between  vorticity  and  stall. 
It  seems  that  vorticity  can  cause  increased  blade  loading. 
Obviously,  if  the  increased  blade  loading  is  of  a  large 
magnitude  and  has  a  sufficiently  long  duration,  this  vorticity 
could  be  the  cause  of  compressor  stall. 

If  it  is  not  possible  to  find  a  logical  correlation  be- 
tween instantaneous  vorticity  and  stall,  it  still  might  be 
possible  to  find  an  empirical  correlation  between  vorticity 
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and  stall  more  satisfying  than  the  relationship  that  has  been 
discovered  between  instantaneous  pressure  and  stall.   This 
would  seem  to  be  a  logical  expectation  since  vorticity  con- 
tains terms  of  a  higher  order  than  pressure  and  may  be  a 
truer  measure  of  distortion. 

These  two  approaches  to  the  interpretation  of  the 
vorticity  maps  both  have  merit,  and  both  should  be  considered 
in  this  paper. 

The  discussion  of  this  topic  will  contain  three  more 
considerations.   First,  it  should  contain  some  type  of 
quantitative  evaluation  of  the  utility  of  the  instantaneous 
vorticity  distribution  in  the  analysis  of  inlet  distortion. 
Secondly  if  the  vorticity  map  proves  to  be  a  useful  tool, 
suggestions  as  to  the  most  beneficial  employment  of  the  map 
in  the  further  study  of  the  problem  will  be  made,  and  some 
work  will  be  done  on  the  development  of  an  inlet  distortion 
factor.   Finally  some  insight  will  be  provided  into  possible 
alternative  paths  to  the  solution  of  the  problem. 
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III.   INTERPRETATION  OF  VORTICITY  MAPS 

A.   EMPIRICAL  APPROACH 

A  computer  program  which  produces  maps  of  instantaneous 
vorticity  at  the  compressor  face  lends  itself  quite  readily 
to  the  empirical  approach  to  the  interpretation  of  these 
maps.   It  produces  a  large  number  of  maps  that  are  easily 
examined.   If  a  pattern  of  vorticity  is  present  in  several 
of  the  events  just  prior  to  stall,  and  not  present  at  other 
times,  it  is  likely  that  the  pattern  represents  the  distortion 
which  was  the  cause  of  compressor  surge.   In  the  examination 
of  these  maps  it  is  necessary  to  consider  the  duration  of  the 
pattern  which  seems  to  be  causing  surge.   If  the  pattern  does 
not  persist  for  a  sufficient  length  of  time,  its  effects  on 
the  compressor  are  likely  to  be  minimized  by  dynamic  suppres- 
sion of  stall.   It  is  possible  to  compute  a  number,  known  as 
the  reduced  frequency,  which  will  provide  some  insight  into 
the  influence  of  dynamic  suppression  on  a  distortion  pattern's 
ability  to  cause  compressor  surge.   Equation  (3),  from  Ref. 
14, 


Reduced  Frequency 


7T  f C     C  =  chord  (3) 


M     (J  =  tip  speed 


is  the  equation  for  reduced  frequency.   For  the  J-S5  the 
length  of  time  that  a  pattern  must  exist  in  order  to  be 
relatively  free  from  the  influence  of  dynamic  suppression  of 
stall  is  approximately  one  millisecond.   This  number,  howevi 
is  more  of  an  order  of  magnitude  approximation  than  a  hard 
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and  fast  rule.   It  could  be  that  a  particularly  strong 
distortion  pattern  with  a  duration  of  slightly  less  than 
one  millisecond  would  be  quite  capable  of  inducing  stall  in 
the  compressor.   For  a  more  complete  explanation  of  dynamic 
suppression  of  stall,  consult  Carta  [15],  Erikson  and  Redding 
[16]  and  Carta  [17] . 

Examination  of  the  maps  generated  by  the  computer  program 
leads  to  the  conclusion  that  the  radial  vorticity  pattern  is 
unlikely  to  yield  an  empirical  correlation  to  stall  or  least 
not  one  that  is  evident  through  visual  inspection.   The 
pattern  of  radial  vorticity  is  generally  quite  random,  and 
the  same  pattern  does  not  persist  for  more  than  a  fraction 
of  a  millisecond.   Figure  6  contains  one  of  the  radial 
vorticity  maps  drawn  by  the  program.   Figure  8  contains  the 
scale  to  be  used  in  examination  of  the  maps.   The  numerical 
values  in  Fig.  7  are  the  values  of  the  constant  normalized 
vorticity  lines  that  are  represented  by  the  boundary  between 
the  two  symbols  on  either  side  of  the  number.   This  scale  is 
the  same  for  all  maps  of  vorticity  that  are  presented  in 
this  thesis. 

The  patterns  of  instantaneous  circumferential  vorticity 
prove  much  more  enlightening.   The  maps  in  Figs.  8-11  are 
patterns  that  appear  prior  to  stall  in  the  first  four  of  the 
eleven  stall  events. 

The  ring  of  large  positive  vorticity  at  the  tip  of  the 
blades  and  extending  nearly  an  entire  360  degrees  seems  to 
be  quite  likely  to  have  caused  stall .   It  is  not  feasible  to 
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FIG.  6  IIlSTAilTMEOUS  RADIAL  VORTICITY  MAP 
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FIG.   8  CASE  I  CIRCUMFERENTIAL  VORTICITY 
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FIG.  9  CASE  II  CIRCUMFERENTIAL  VORTICITY 
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FIG.  11  CASE  IV  CIRCUMFERENTIAL  VORTICITY 
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examine  every  map  in  a  stall  event,  due  to  the  amount  of 
computer  time  required  for  the  production  of  these  maps; 
each  set  of  three  maps  requires  50  seconds  of  computer  time. 
With  2048  sets  of  three  maps  per  stall  event,  this  would 
amount  to  2  8  hours  of  computer  time  for  each  event.   A 
random  sampling  of  vorticity  maps  was  taken,  and  the  ring 
pattern  does  not  appear  except  at  such  a  time  that  it  could 
be  the  cause  of  stall.   The  similarities  between  these  maps 
are  striking.   A  distortion  index  based  on  this  pattern 
could  yield  some  insight  into  how  the  compressor  stalls. 
It  would  appear  from  this  pattern  that  the  tips  of  the  blades 
are  the  cirtical  point  of  the  compressor  for  the  J-85. 

B.   ANALYTICAL  APPROACH  [Ref.  8] 

There  has  been  a  large  amount  of  work  done  for  predicting 
the  effects  of  nonuniform  flows  on  the  performance  of  a 
cascade.   Reference  8  contains  an  excellent  summary  of  some 
of  the  results  of  this  work  as  well  as  sample  calculations 
of  the  effects  of  secondary  flow  on  blade  loading.   Further 
insight  into  these  effects  can  be  gained  by  consulting  Ref.  14. 

The  primary  effects  of  radial  vorticity  can  be  understood 
easily.   Reference  18  states  that  an  airfoil  can  be  represened 
by  a  reentrant  vortex.   It  is  fairly  easy  to  see  that  the 
blades  in  a  compressor  could  be  approximated  by  a  group  of 
negative  radial  vortices.   (See  Fig.  12).   It  is  then  a  logical 
supposition  that  positive  radial  vorticity  would  tend  to  de- 
crease the  blade  loading  and  negative  radial  vorticity  would 
tend  to  increase  blade  loading,  as  in  Fig.  13.   The  trouble 
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with  this  observation  is  that  radial  vorticity  in  a  rotor 
is  an  unsteady  phenomenon.   Table  II  classified  types  of 
vorticity  in  a  rotor  and  stator  as  steady  or  unsteady. 


FIG.  12  EQUIVALENCE  OF  BLADE  CIRCULATION  AMD  VORTICITY 


FIG.  13  ADDITIONAL  LIFT  DUE  TO  RADIAL  VORTICITY 


T/\BLE  ii.  POSSIBLE  COMBINATIONS  IN  A  ROTOR  [Ref.  8] 
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The  influence  of  circumferential  vorticity  on  blade 
loading  is  not  so  clear.   In  Ref.  8  the  secondary  effects 
of  one  vorticity  pattern  were  calculated  with  the  following 
results.   It  was  found  that  circumferential  vorticity  could 
cause  additional  blade  loading  through  its  secondary  flow 
effects.   The  vorticity  pattern  for  which  these  effects  were 
calculated  is  remarkably  similar  to  the  pattern  that  the 
empirical  correlation  has  tentatively  singled  out  as  the 
cause  of  compressor  surge.   This  pattern  is  depicted  in 
Fig.  14.   In  Fig.  15  the  induced  downwash  from  this  vorticity 
pattern  is  pictorialized.   If  the  blade  tips  were  highly 
loaded,  the  additional  lift  caused  by  the  circumferential 
vorticity  definitely  would  be  large  enough  to  initiate  stall. 
It  seems,  therefore,  that  there  is  a  logical  relationship  be- 
tween a  ring  of  large  positive  vorticity  at  the  tip  and  stall 
that  concurs  with  the  purely  empirical  correlation  which  was 
pointed  out  in  the  previous  section.   For  a  thorough  expla- 
nation of  the  method  of  calculation  of  these  secondary  effects, 
consult  Ref.  8  or  Ref.  1. 


FIG,  I'!  VORTICITY  PATTER! 
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FIG,  15  INDUCED  DOWflWASH   [REF.  8] 


FIG,  16  REG  I OI1S  OF  COflPRESSOR  FACE 
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IV.   A  DISTORTION  INDEX 

A.   DEVELOPMENT  OF  THE  INDEX 

From  the  two  approaches  to  the  interpretation,  it  has 
been  seen  that  stall  in  the  J-85  is  probably  caused  by  a 
type  of  distortion  which  is  characterized  by  a  large  nega- 
tive vorticity  at  the  midspan  and  a  large  negative  vorticity 
at  the  tip.   Both  interpretations  point  to  this  conclusion, 
making  it  fairly  certain  that  this  path  of  investigation  will 
bear  fruit.   In  order  to  investigate  the  problem  further  and 
to  determine  if  the  conclusion  that  stall  is  caused  by  that 
one  particular  vorticity  pattern  is  correct,  it  would  be 
helpful  to  develop  a  distortion  index  based  on  the  assumption 
that  that  pattern  causes  stall. 

The  distortion  index  will  have  certain  advantages  over 
the  maps.   The  most  important  of  these  is  that  it  will  take 
less  computer  time  and  will  facilitate  the  processing  of  more 
data.   This  will  make  it  possible  to  evaluate  data  for  which 
no  map  was  produced  and  to  check  the  conclusion  made  in  the 
interpretation  of  the  maps. 

The  index  should  fulfill  certain  requirements.   It  should 
be  large  when  the  vorticity  pattern,  tentatively  labeled  as 
the  stall  producing  pattern,  is  present  and  smaller  when  it  is 
not.   It  should  be  large  only  if  the  pattern  persists  for  one 
millisecond  or  more,  since  if  it  does  not  persist  for  that 
long,  it  is  unlikely  that  it  would  have  produced  stall. 
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Since  this  is  the  first  index  ever  devised  based  on  vorticity, 
there  is  no  guideline  as  to  what  is  large  and  what  is  small. 
This  is  an  advantage  in  that  it  does  not  matter  how  large  the 
index  is.   It  is  only  the  relative  size  at  different  points 
in  time  that  is  of  interest.   However,  it  will  be  possible  to 
get  some  feel  for  the  size  of  this  index  which  causes  stall 
through  the  comparison  of  the  index  for  different  stall  events, 

The  development  of  the  index  will  now  be  described. 
Through  inspection  it  was  seen  that  the  positive  ring  of 
vorticity  was  generally  at  r '  of  greater  than  0.75,  and  the 
negative  ring  of  vorticity  was  at  r'  of  less  than  0.75.   The 
boundary  layer  had  dominated  the  vorticity  at  r1  greater  than 
0.95.   At  r'  less  than  0.5,  there  was  very  little  regularity 
to  the  pattern.   In  order  to  reduce  computer  time,  it  was 
decided  to  employ  the  observations  above  and  to  make  the  in- 
dex a  simple  arithmetic  summation  with  as  little  testing  for 
magnitude  as  possible.   To  this  end,  the  vorticity  map  was 
separated  into  three  sections  as  shown  in  Fig.  16.   Region  I 
at  r'  less  than  0.5  was  ignored.   In  region  II,  the  vorticity 
at  radii  of  0.55,  0.60,  0.65,  and  0.70  was  found  for  0  of  0 
to  360  at  steps  of  0.5°,  and  summed  yielding  equation  (4). 

___        7/7 

If  the  vorticity  at  any  point  in  this  region  was  greater  than 

0,  it  was  rejected.   In  region  III,  the  positive  vorticity  at 

radii  of  0.8,  0.85,  0.90,  and  0.95  was  foun<i  and  summed.   In 

this  case  negative  numbers  were  rejected.   This  yields 

equation  (5) . 
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U5 
Facta (t)=  /L  fyo*T \!ioitiii.5)-x&(ctt.9s;NK.s)  t   vfc«T  (•'<>>*'*•*')  +^(.1^^)]  (5) 

The  distortion  factor  is  then  the  difference  between  these 
two  sums,  (equation  6) . 

FftCT  (t)  -  facT^X  Cc)  -  f7»cTJ  (t)  (6) 

At  this  point,  a  distortion  factor  was  known  at  each 
instant  of  time.   If  the  ring  of  positive  vorticity  at  the 
tip  is  present  along  with  the  ring  of  negative  vorticity  at 
the  hub,  this  index  should  be  a  large  number.   However,  this 
number  carries  no  information  about  the  duration  of  the 
pattern.   In  this  simple  development  the  time  factor  was  in- 
corporated in  the  following  manner.   The  distortion  factor 
for  each  case  was  multiplied  by  the  distortion  factor  0.625  ms, 
before  it  and  0.6  25  ms.  after  it.   In  this  manner  a  new  dis- 
tortion factor  is  calculated  that  reflects  the  vorticity  pat- 
tern as  seen  through  a  "window"  that  is  1.2  5  milliseconds 
long.   This  is  a  simple  matter  to  program  since  the  cases  used 
are  all  0.625  milliseconds  apart.   The  distortion  factor 
becomes 


K  Dt*.T  ~   Fact  (.  t  -  ..** S"  *n t)  x  \-,\ cr  t±)    x   h^cr(t  t.ws}t 7 

In  this  manner  a  particular  distortion  must  last  at  least 
1.0  millisecond  in  order  to  yield  a  large  value  of  the 
distortion  factor. 
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B.   TEST  OF  THE  DISTORTION  FACTOR 

The  distortion  factor  was  inserted  into  the  computer 
program  and  tested  for  several  stall  events.   A  graphical 
representation  of  the  results  of  each  of  the  first  four 
events  is  presented  in  Figs.  17-20. 

It  can  be  seen  that  the  distortion  factor  is  quite 
effective  for  three  of  the  four  cases.   In  the  fourth  case, 
however,  it  gives  no  indication  whatsoever  of  stall.   This 
could  be  due  to  one  or  all  of  several  influences.   The  index 
fails  to  account  for  several  things.   It  could  be  that  the 
index  becomes  large  at  the  wrong  time  due  to  the  fact  that 
there  may  be  a  ring  of  positive  vorticity  at  the  tip  and 
negative  vorticity  at  the  midspan  without  the  increase  in 
lift  if  they  do  not  coincide  as  shown  in  Fig  21.   If  this 
is  the  case,  they  could  cause  the  index  to  be  large  without 
inducing  stall.   Also,  there  may  be  other  influences,  that 
are  at  the  present  obscure,  which  play  a  part  in  the  stall  of 
a  compressor.   Finally,  it  could  be  that  a  mild  distortion 
can  cause  stall  if  it  persists  for  long  enough. 

Another  interesting  point  about  the  test  is  that  for  two 
of  the  three  successful  tests  the  peak  of  the  distortion  index 
which  seems  to  have  initiated  stall  is  of  the  same  order  of 
magnitude,  while  the  third  is  much  smaller.   This  seeming 
discrepancy  could  provide  some  insight  into  the  mechanism  of 
stall.   The  operating  points  for  the  four  cases  were  plotted 
on  a  performance  map  for  the  J-0  5,  and  the  steady-state 
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operating  point  for  the  case  in  which  a  smaller  peak  in 
distortion  index  initiated  stall  is  much  closer  to  the  surge 
line  than  the  other  two  successful  tests.   This  may  be  a 
significant  observation. 

In  any  case,  if  the  conclusions  that  lead  to  the  dis- 
tortion index  are  correct,  then  the  initial  level  of  dis- 
tortion is  not  the  same  for  all  operating  conditions  of  the 
compressor.   Also,  there  is  something  about  that  particular 
operating  point  for  the  fourth  case  that  makes  the  compressor 
more  susceptible  to  stall. 


FIG,  21  VORTICITY  PATTERN  PRDUCES  LARGE  INDEX, 
HO  STALL 
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FIG,  -28  DISTORT  I  Oil  INDEX,  CASE  II 
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FIG,  10  DISTORTION  1,  DEX,  CASE  III 
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FIG,  20  DISTORTION   1  I       ,  CASF.  IV 
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V.   SUMMARY  AND  REVIEW 

In  review  it  is  necessary  to  put  the  findings  of  this 
thesis  in  perspective.   In  order  to  test  the  usefulness  of 
the  instantaneous  vorticity  approach  to  the  problem  of  in- 
let distortion  in  axial  flow  compressors,  the  pressure 
distribution  at  the  compressor  face  must  be  converted  to  a 
vorticity  distribution  at  the  compressor  face.   A  computer 
program  was  developed  which  performed  this  task  with  reason- 
able accuracy. 

Once  the  program  was  operational  it  was  used  to  convert 
instantaneous  pressure  data  for  the  J-85-GE-13,  which  was 
taken  in  the  10  x  10  foot  supersonic  wind  tunnel  at  NASA 
Lewis  Research  Center  in  Cleveland,  Ohio,  into  vorticity 
maps  of  radial  and  circumferential  vorticity  at  the  compres- 
sor face.   The  map  of  radial  vorticity  offered  little  infor- 
mation, but  a  marked  empirical  relationship  was  found  between 
a  circumferential  vorticity  pattern  and  stall.   A  crude 
instantaneous  distortion  index  was  defined  using  this  empirical 
relationship  as  a  criterion.   The  results  of  these  tests  show 
that  there  is  indeed  an  empirical  relationship  between  instan- 
taneous vorticity  and  stall. 

The  tests  of  the  distortion  factor  also  showed  that  this 
pattern  of  vorticity  is  not  the  only  mechanism  that  is  capable 
of  causing  stall. 
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VI.   CONCLUSIONS 

The  primary  purpose  was  to  establish  the  validity  of 
the  vorticity  approach  to  the  problem  of  instantaneous  dis- 
tortion in  a  jet  inlet.   This  approach  to  the  problem  has 
proved  to  be  quite  useful  in  understanding  the  effects  of 
distortion  on  the  stability  of  a  compressor.   The  computer 
program  which  converts  the  pressure  distribution  at  the 
compressor  face  to  the  vorticity  pattern  at  the  compressor 
face  will  be  instrumental  in  further  research  along  this 
line  of  reasoning. 

While  the  crude  distortion  factor  developed  in  this  work 
is  fairly  successful  in  the  prediction  of  stall,  and  could 
most  likely  be  refined  to  be  even  more  effective,  it  is  not 
the  purpose  of  this  line  of  research  to  determine  an  empiri- 
cal instantaneous  distortion  pattern  for  the  J-85.   The  data 
for  the  J-85  were  used  to  show  that  the  vorticity  approach 
to  the  problem  is  valid.   The  overall  purpose  of  the  re- 
search is  to  develop  a  universal  inlet  distortion  factor 
that  could  be  employed  for  a  wide  variety  of  engines  and 
airframes.   Further  research  should  be  along  these  lines. 

It  may  be  possible  to  make  progress  toward  a  universal 
inlet  distortion  factor  through  the  calculation  of  the 
theoretical  secondary  flow  effects  of  vorticity  as  it 
passes  through  the  compressor.   In  further  work  the  radial 
and  axial  vorticity  must  be  incorporated  into  the  analysis. 
Other  factors  that  must  be  considered  in  the  analysis  include 
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blade  twist  and  type  of  blading.   Finally,  as  stated  in 
Ref.  2,  other  types  of  distortion,  such  as  temperature 
distortion,  must  be  considered.   The  vorticity  approach  to 
the  problem  of  inlet  distortion  is  very  promising  and  may 
provide  the  reliable  universal  distortion  criteria  that 
industry  lacks  today. 
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APPENDIX  A 
EXPLANATION  OF  DATA 

The  data  used  by  the  computer  program  are  the  same  data 
which  were  used  in  Ref.  5.   They  were  taken  in  the  10  x  10 
foot  supersonic  wind  tunnel  at  NASA  Lewis  Research  Center  in 
Cleveland,   They  consist  of  eleven  different  stall  events. 
For  each  stall  event  the  250  milliseconds  of  analog  data 
prior  to  stall  were  digitized  at  8  0  00  points/second/channel. 
A  schematic  of  the  digitizing  system  is  in  Fig.  A-l.   Each 
of  the  digitized  values  consists  of  only  the  fluctuating 
component  of  pressure,  so  the  steady-state  values  of  pressure 
just  prior  to  stall  must  be  used  in  order  to  find  the  absolute 
total  pressure  at  each  point.   The  steady-state  value  was  not 
used  alone,  however.   Instead  a  reference  pressure  which  in- 
cludes a  correction  for  the  leakage  of  the  high  pass  filter 
was  added  to  the  fluctuating  component  to  obtain  absolute 
values  of  the  stagnation  pressure.   The  reference  pressure 
was  determined  in  the  following  manner.   Analog  data  for  all 
thirty  channels  were  filtered  at  50  hertz.   A  drawing  of  the 
compressor  face  with  the  positions  of  all  the  probes  is  shown 
in  Fig.  A-2.   A  five  second  section  of  data  prior  to  stall 
was  digitized  at  200  points/second/channel.   The  thousand 
values  for  each  channel  were  then  averaged  to  find  delta  p 
average  for  each  channel,  where  delta  p  average  is  equal  to 
the  filter  leakage.   The  values  of  Pref  and  dP/dt  then  were 
calculated  using  equations  A-.l  and  A-2. 
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FIG,  A-2  PLACEMENT  OF  PROBES 
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FfitF  lsri5Aiiy-iTAT£       -       A   I  /W  6  (A    1) 

Pao6     ~    Uef     f     l~LucraftTif06      C6MPotdeeJT  (A-2) 

For  each  event  the  time  of  engine  stall  was  known  as  well 
as  the  time  of  earliest  indication  of  stall  inside  the 
compressor.   This  time  of  earliest  indication  of  stall  inside 
the  compressor  is  the  time  at  which  the  interstage  dynamics 
indicated  that  the  engine  was  going  to  stall  and  is  accurate 
only  to  the  nearest  millisecond.   A  synopsis  of  these  times 
for  each  stall  event  is  in  Table  A-I.   Also,  a  list  of 
steady-state  operating  conditions  for  each  stall  event  is 
contained  in  Table  A-II. 

If  a  more  complete  explanation  of  the  instrumentation 
used  is  desired,  consult  Ref.  5. 
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APPENDIX  B 
ACCURACY  OF  THE  COMPUTER  PROGRAM 

When  a  computer  program  of  this  complexity  is  used,  it 
is  mandatory  that  the  accuracy  of  the ' program  be  verified. 
In  the  calculation  of  the  vorticity  distribution  at  the  com- 
pressor face  there  are  six  quantities  involved:   r ' ,  $  , 
normalized  axial  velocity,  normalized  total  pressure,  and  the 
derivatives  of  the  normalized  pressure  in  the  r"  and  theta 
directions.   The  distance  r'  is  actually  one  of  the  defining 
values  for  any  point  on  the  compressor  face,  and  gamma  is 
merely  the  ratio  of  specific  heats  which  is  a  well  known  func- 
tion of  temperature.   The  four  remaining  quantities  are  cal- 
culated by  the  program  as  a  function  of  r'  and  0  that  is  pro- 
vided by  the  interpolation  scheme.   Since  the  normalized 
axial  velocity  is  a  function  of  two  measured  quantities,  the 
average  total  pressure  and  the  average  static  pressure,  and 
the  total  pressure,  it  is  obvious  that  if  the  interpolation 
scheme  accurately  computes  the  total  pressure  at  each  point, 
then  the  normalized  axial  velocity  also  will  be  correct. 
Thus  we  sec  that  the  accuracy  of  the  entire  program  hinges 
on  the  accuracy  of  the  interpolation  scheme  which  calculates 
P '  01   ,  and  C'S    .   The  interpolation  scheme  is  a  cubic  spline 
fit  in  two  variables .   It  fits  a  cubic  polynomial  to  each 
three  points  in  one  dimension  and  then  matches  the  slope  of 
the  curve  at  the  points  that  are  common  to  both  polynomials. 
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There  are  two  ways  to  test  the  accuracy  of  this  inter- 
polation scheme.   One  way  is  to  set  up  a  dummy  function  of  r1 
and  Q,    and  produce  a  grid  of  point  values  of  this  dummy 
function.   Then,  using  the  interpolation  scheme,  predict  the 
value  of  the  dummy  function  and  its  two  derivatives  at  other 
points  in  between  the  grid  values  and  compare  these  values 
to  the  true  values  of  the  dummy  function  which  are  known 
since  the  dummy  function  is  known  explicitly.   The  second  way 
is  to  plot  the  known  values  of  the  pressure  on  a  graph,  con- 
nect the  points  with  a  curve,  and  see  if  the  interpolated 
values  fall  on  the  curve.   The  slope  of  the  curve  can  be 
compared  with  the  partial  derivatives  that  were  calculated. 
The  results  of  the  first  method  of  testing  are  in  Table  B-l. 
In  Figs.  B-l  -  B-5  the  results  of  the  second  method  are 
presented. 

The  interpolation  scheme  was  found  to  have  a  maximum 
error  of  one  percent  for  the  value  of  pressure  in  this  test, 
while  the  values  of  the  partial  derivatives  were  in  error  by 
a  maximum  of  four  percent. 

Once  the  accuracy  of  the  interpolation  scheme  has  been 
verified,  the  only  remaining  possibility  of  error  is  the  sub- 
routines that  print  the  map  of  the  vorticity  distribution. 
In  Table  B-II  the  values  of  vorticity  calculated  on  a  desk 
calculator  are  compared  with  the  symbols  that  were  printed 
at  these  points  by  the  computer  while  drawing  the  maps. 
The  positions  of  these  points  are  indicated  on  the  vorticity 
map  in  Fig.  B-6. 


56 


TABLE  B-I 
PREDICTED  VALUES  OF  A  FUNCTION  VS  ACTUAL  VALUES 


Actual 

1918.2 

1927.3 

1947.5 

1982.0 

2033.9 

2106. 3 

1723.9 

1693.9 

1657.9 

1615.9 

1567.9 

1513.9 

1382.6 

1313.5 

1221.3 

1102.8 

954.9 

774.5 

941.8 

872.7 

780.5 

661.9 

514.1 


Predicted   Actual 
1920.5      -  97.49 


Predicted   Actual 
-100.3        45.66 


1929.4 

1950. 0 

1984.9 

2036.8 

2111.3 

1723.0 

1693.3 

1657.1 

1614.9 

1566.8 

1512.1 

1382.1 

1314.2 

1222.0 

1103. 8 

956.9 

774.9 

941.3 

873.4 

781.2 

662.9 

516.0 


-120.0 
-152.5 
-196.7 
-254.3 
-327.3 
-265.1 
-310.3 
-375.2 
-463.6 
-578.9 
-724.9 
-378.1 
-400.7 
-433.1 
-477.3 
-535.0 
-608.0 
-463.7 
-441.2 
-408.7 
--364.5 
-306.9 


-122.9 
-156.1 
-201.1 
-259.6 
-355.0 
-265.2 
-309.6 
-374.8 
-463.2 
-578.1 
-726.5 
-378.7 
-401.0 
-4  3  3.7 
-47  8.1 
-535.7 
-610.1 
-4  6  3.1 
-440.9 
-408.2 
-363.8 
-306.1 


141.5 

268.6 

426.8 

616.2 

836.8 

-269.9 

-329.9 

-389.9 

-449.9 

-509.9 

-569.9 

-585.6 

-801.5 

-1048.6 

-1326.8 

-1636.2 

-1976.8 

-585.6 

-8  01.5 

-1048.6 

-1326.8 

-1636.2 


Predicted 

43.59 

144.6 

271.0 

431.7 

617.6 

856.6 

-267.5 

-331.9 

-391.1 

-452.2 

-511.6 

-576.3 

-574.4 

-802.1 

-1044.5 

-1324.4 

-1625.8 

-1990.0 

-574.4 

-802.1 

-1044.5 

-1324.4 

-1625.8 
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FIG.  B-2     PRESSURE  vs  TIIETA,  r'=0.59 
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TABLE  B-II 
TEST  OF  PRINT  ROUTINE  SYMBOLS 


r. 

0 

W 
Q 

Symbol 
Printed 

Symb 
Shou 

ol  Which 

Id  be  Printed 

.41 

20° 

1.055 

- 

- 

.59 

60° 

-6.15 

H 

H 

.73 

100° 

-1.  09 

E 

E 

.59 

140° 

.59954 

* 

* 

.41 

180° 

-1.64 

H 

H 

.85 

180° 

-1.89 

H 

H 

.95 

180° 

-.1986 

(D 

CD 

.73 

180° 

-2.1 

H 

H 

.85 

340° 

-.296 

CD 

CD 

.95 

340° 

.1657 

0 

O 
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